Real-time observation of discrete Andreev tunneling events 
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We provide a direct proof of two-electron Andreev transitions in a superconductor - normal metal 
tunnel junction by detecting them in a real-time electron counting experiment. Our results are 
consistent with ballistic Andreev transport with an order of magnitude higher rate than expected 
for a uniform barrier, suggesting that only part of the interface is effectively contributing to the 
transport. These findings are quantitatively supported by our direct current measurements in single- 
electron transistors with similar tunnel barriers. 
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Electronic transport across a boundary between con- 
ductors with dissimilar carriers is a non-trivial process. 
Of particular interest in this respect is the transport 
through a superconductor - normal metal interface that 
at low energies is dominated by Andreev reflection 
[lol | , where a Cooper pair in a superconductor is converted 
into two electrons in the normal metal or vice versa. Here 



we employ electron counting techniques llMl7j| to detect 
these Andreev events. Since the observed rate depends 
on the coherence of the two electrons involved in the tran- 
sition, we obtain, as a result, a fingerprint of the junction 
electrodes and the tunnel barrier. 

The techniques used for observing individual electrons 
are based on the Coulomb blockade effect where the elec- 
trostatic energy of a small metallic island changes no- 
ticeably when only one elementary charge e is placed 
on or removed from it. In the present experiment, we 
employ an isolated single-electron box where a supercon- 



ducting island is connected to a normal metal one 17|, 
but neither of these two is connected galvanically to the 
external circuitry. The electron tunneling rates between 
the islands are then sufficiently low to be monitored by 
low-frequency electrometry and are described in detail by 
relatively simple theoretical considerations 1^ isl. We 
use a single-electron transistor (SET) [uHiillllll-iil 
as an ultrasensitive electrometer. With charge sensitiv- 
ity as good as 10~^e/'\/Hz 22-24], it is capable to detect 
individual electrons with high precision. In Fig. [U we 
show a micrograph of our sample fabricated by standard 
e-beam processing. 

The tunneling rates and resulting charge distribution 
between the two islands of the isolated box can be ad- 
justed with an offset charge induced by a gate voltage. 
The electrostatic energy of a state with n excess elec- 
trons on one of the islands is given by = Ec{n — rig)'^ , 
where is the charging energy for individual electrons, 
and Ug is the normalized offset charge that can be viewed 
as the polarization charge on the gate capacitor and de- 
termines the energetically preferred way to occupy the 
different charge states n (Ref. 20]). In Fig. [2] (a), the 
two extreme cases are shown. In the Coulomb blockade 




FIG. 1. (color online) Scanning electron micrograph of the 
measured structure and the schematic layout of the measure- 
ment set-up. The isolated electron box consists of two metal- 
lic islands, seen as 25 fim long rectangles (colored red). They 
are connected to each other by a normal metal - insulator - 
superconductor (NIS) tunnel junction. Tunneling of electrons 
through the junction is monitored with a DC SET electrom- 
eter (in blue) coupled capacitively to one of the box islands. 
It is voltage biased, and the gate voltage of the SET sets the 
electrometer to a charge sensitive point, and the gate volt- 
age of the box tunes the energy levels of the charge states in 
it. The NIS junction (top) and detector (bottom) are shown 
magnified on the left side of the main micrograph. 



regime for single electrons, Ug is an integer and the state 
n — Ug has the minimal energy i?„ = 0. To enter an 
excited state, one electron can either tunnel into or out 
of the island (dotted black lines with arrows in Fig. [2] 
(a)) but energy Ec has to be provided for the tunneling 
electron in addition to the Cooper pair breaking energy 
equal to or larger than the superconducting energy gap 
A [25]. In the other extreme, at degeneracy with half- 
integer Ug, two electron states differing by charge e have 
equal minimal energy and hence are equally populated. 
The tunneling rate between them is higher than in the 
Coulomb blockade regime as no extra energy for charg- 
ing is needed. For Andreev reflection (solid red lines with 
arrows in Fig. [5] (a)), the energy cost of charging is cal- 
culated similarly, but now the initial and flnal states are 
separated by two electrons, and the energy cost of break- 
ing a Cooper pair is avoided since complete pairs tunnel 
at once. 

In the experiment (see supplemental material for ad- 
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FIG. 2. (color online) Energy levels of the various charge states and typical observed time traces of the current through the 
detector, (a) Low lying levels of the box in the Coulomb blockade (upper panel) and at degeneracy (lower panel). Dotted 
(black) and solid (red) arrows indicate one- and two-electron processes, respectively, (b) Measured time traces of the detector 
current showing the charge state of the electron box as a function of time. The leftmost panel presents the case of Coulomb 
blockade. The rightmost panel depicts the opposite limit where the two charge states are equal in energy (degeneracy). The 
trace in the centre is taken halfway between these two cases, (c) Grey sections of the traces of (b) zoomed. Vertical arrows 
indicate two-electron events. 



ditional details [26j), we measured time traces of the de- 
tector current at various biasing conditions of the box. 
With all the other parameters fixed, the detector gate 
voltage was adjusted to maximize charge sensitivity and 
dynamic range. The observed current jumps (Fig. 2b 
and c) are attributed to the tunneling events between 
the two islands. The switching rate depends on the gate 
voltage of the box and hence on its charge state, being 
lowest in the Coulomb blockade regime (leftmost panel) 
and highest at degeneracy (rightmost panel). The events 
observed in the traces (grey regions of Fig. [2] (b) zoomed 
in [2] (c)), indicate that individual electrons tunnel be- 
tween the islands: In the Coulomb blockade regime they 
hop rarely from the lowest (n — 0) to the higher (n = ±1) 
energy states and back, while at degeneracy the electrons 
tunnel frequently between the two lowest states and only 
occasionally the system enters a higher lying level. More 
interestingly, the traces show also the coincident events, 
pointed out by the vertical arrows, where two electrons 
appear to tunnel simultaneously. In the following, we 
show that most of these events represent Andreev tun- 
neling. 

Due to the finite measurement bandwidth, limited to 
1 kHz by the dc readout of the electrometer, events re- 
sembling the two-electron Andreev tunneling could in 
principle arise from almost coincidental tunneling of two 
independent quasiparticles. To assess this option, we 
recorded time traces for several minutes at each gate off- 
set value. From the traces we determined the distribution 
of the time t spent in the state n = before a transition 
took place. In Fig. [3] (a) we show such a distribution on 



the left for Coulomb blockade {ug = 0), in the centre for 
Ug — 0.25, and on the right near degeneracy {ug — 0.45). 
A direct transition (Andreev tunneling) between states 
n = ±1 contributes here as essentially a,t = event since 
the time separation between the two electrons tunneling 
in the Andreev process should be on the order of h/A, 
which is many orders of magnitude smaller than the time 
scales relevant in Fig [3] (a). Overall, the distribution is 
exponential as we have Poisson distributed one-electron 
tunneling processes. However, at small lifetimes in the 
Coulomb blockade regime, the data point indicated by 
the horizontal arrow does not follow the exponential de- 
pendence and corresponds to excessively many events. 
This clear separation of the short-lifetime events from 
the one-electron transitions shows that the majority of 
these events are not coincidental one-electron tunneling 
but rather two electrons tunneling concurrently. When 
the box offset charge is adjusted closer to degeneracy, 
the anomalous datapoint gradually merges to the rest of 
the data, in accordance with its interpretation in terms 
of Andreev transitions, since the energy cost of the two- 
electron tunneling increases. Eventually at degeneracy, 
almost no two-electron events are observed. We empha- 
size that the charging energy should be small, Ec < A, 
for Andreev tunneling to occur, since for large Ec, it is 
not energetically favourable [17|. 

For quantitative analysis we counted the number of 
one- and two-electron events, and normalized those by 
the time spent in the initial state to obtain the corre- 
sponding tunneling rate. The results for both processes 
at the temperature of 60 mK are shown in Fig.[3] (b). The 
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FIG. 3. (color online) Lifetime distributions and tunneling 
rates, (a) Lifetime distributions of the n = charge state 
in the Coulomb blockade, at an intermediate gate position 
and close to degeneracy, from left to right. The correspond- 
ing gate positions are indicated by vertical lines in (b). The 
distribution is exponential as expected for a Poisson process. 
However, at small lifetimes an anomaly, pointed out with the 
horizontal arrows, is observed as two-electron events are be- 
ing counted as fast events. The anomaly vanishes as one ap- 
proaches degeneracy since the energy cost is high for the two- 
electron events, (b) One-electron (open symbols) and two- 
electron (solid symbols) tunneling rates estimated from the 
counted events at the base temperature. Red circles and blue 
triangles denote the tunneling rates out of states n = — 1 and 
1. Dark blue diamonds and dark yellow squares denote tun- 
neling rates out of state to states —1 and 1. Dark green 
solid lines present calculated Andreev rates with C = 4 • 10~^. 



parameter values which determine the first order tunnel- 
ing rates in the box junction, the tunneling resistance 
Rt = 2 Mil, superconducting energy gap A = 216 /^eV, 
and the charging energy Ec = 0.2 A, were determined 
from the one-electron tunneling rates measured at tem- 
peratures ranging from 60 mK to 200 mK as in Ref. ■ 

To calculate the two-electron rates, we assume ballistic 
model of the tunnel barrier, which fits well our experi- 
ments with small junctions. Apart from the parameter 
values obtained from one-electron tunneling, the only pa- 
rameter to be determined for the two-electron tunneling 
is its overall magnitude controlled by C = (?i/^Te^)/A/', 
where M is the effective number of the conduction chan- 
nels in the junction, M = A/Ach [llj- The junction area 
A = 40 nm x 35 nm is estimated from the scanning elec- 
tron micrograph of the tunnel junction, and the area of 
one conduction channel Ach = 30 nm^ is determined in- 
dependently from measurements of the SET current de- 
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FIG. 4. (color online) Detection of Andreev current in SETs. 
Colored regions show measured subgap current for devices 
with charging energies Ec/ks ~ 2.3, 1.9, 1.5 and 0.86 K 
from top (yellow) to bottom (orange) at all gate offset val- 
ues. The insets depict the larger scale measurements for 
samples with the highest and lowest charging energies and 
a scanning electron micrograph of one of the measured SETs. 
Solid black lines are fitted theoretical curves at degeneracy 
(maximal current) and in Coulomb blockade (minimal cur- 
rent) when Andreev tunneling is taken into account. Dotted 
blue lines present fits excluding Andreev processes. Tunnel- 
ing resistances of the samples were Rt ~ 129, 78, 55 and 
31 kil in the order of decreasing Ec, and superconducting 
gap A = 216 ^eV for all of them. Fitting parameter for the 
magnitude of the Andreev processes was C ^ 4 - 10-*. (Note 
that to increase the measured signal the SET junctions were 
oxidized less to have larger specific conductance than the box 
junction.) Open (solid) circles present the expected thresh- 
olds, eV — ±2Ec, (eV = ±4_Bc) for Andreev tunneling at 
degeneracy (Coulomb blockade). 



scribed in detail below. The results of the calculation of 
the Andreev rates with these parameter values are plot- 
ted as sohd lines in Fig.[3](b), and show good agreement 
with the measured rates. For completely uniform tun- 
nel junction, however, the area of a conduction channel 
should be about Ach — 2 nm^, see supplemental ma- 
terial [26]. We note that the measured value is about 
one order of magnitude higher than the theoretical esti- 
mate, which agrees with the conclusions in, e.g. [jl. [lo|. 
for larger tunnel junctions that exhibit diffusive Andreev 
tunneling. We attribute this to the imperfections of the 
barrier, with only a small portion of the junction area 
dominating the transport. This fact can be understood 
by noting that the standard parameters of the aluminum 



oxide tunnel barriers [26|, |28|, |29j imply that, for exam- 



ple, variation of the barrier thickness by one atomic layer 
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(i.e. by about 0.3 nm) results in the specific conductance 
change by about a factor of 10. This means that even 
relatively small fluctuations of the barrier thickness of 
the magnitude of one to two atomic layers in compari- 
son with the average barrier thickness of 6 atomic layers 
can reduce the effective area of the region dominating the 
barrier transparency to about 10% of the total junction 
area. 

As an additional way of extracting parameters of 
the two-electron tunneling, we measured current-voltage 
characteristics of basic SET structures with supercon- 
ducting leads and a normal metal island. A micrograph 
of one of the devices is shown in the inset of Fig. |4] (a) . 
All four SETs of Fig. |4] were fabricated in the same batch 
so that only the junction area (Ec) varies. Here An- 
dreev tunneling has to be extracted from the measured 
data containing a contribution also from one electron pro- 
cesses. The SET was biased at voltage V and the cur- 
rent / flowing through the device was measured. At each 
bias voltage, the gate offset of the island was varied so 
that the maximal and minimal current, corresponding 
to degeneracy and Coulomb blockade, respectively, were 
observed. The measured current values between these 
extremes are shown as the colored regions in Fig. |4] (b) 
for four devices having different charging energies. For 
the highest charging energy (top), we have no current 
flowing in the subgap region eV < 2 A, apart from the 
thermally activated single-electron tunneling. When the 
charging energy is lowered, the onset of the Andreev cur- 
rent penetrates into the subgap region with the threshold 
at eV — 2Ec, seen both in the experimental data and in 
the flts including Andreev processes (solid black lines) , as 
the energy cost for two electrons tunneling into the island 
is lowered. With the high charging energies, Ec > A, it 
is sufficient to take into account only one-electron tun- 
neling in the fits (dotted blue lines). On the large current 
scale, the simulations with and without Andreev current 
are indistinguishable as shown in the insets for samples 
with the highest and the lowest charging energy. In the 
numerical simulations for all the samples, we obtain as 
the only fit parameter the area of one conduction chan- 
nel to be Ach = 30 nm^. This value was adopted to the 
analysis of the counting data above. 

In conclusion, we have detected in time-domain in- 
dividual Andreev transitions allowing us to distinguish 
higher order tunneling from the usual one-electron pro- 
cesses. These counting experiments are consistent with 
our direct current measurements, and they both indicate 
that the tunnel barriers are non-uniform. 
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